Optimal control of light pulse storage and retrieval 
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We demonstrate experimentally a procedure to obtain the maximum efficiency for the storage and 
retrieval of light pulses in atomic media. The procedure uses time reversal to obtain optimal input 
signal pulse-shapes. Experimental results in warm Rb vapor are in good agreement with theoretical 
predictions and demonstrate a substantial improvement of efficiency. This optimization procedure 
is applicable to a wide range of systems. 
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Mapping of quantum states between light and mat- 
ter is a topic of great current interest pj 0, Q. One 
of the leading approaches to realizing this capability 
is the storage of light in ensembles of radiators (warm 
atoms, cold atoms, impurities in solids, etc.) using a dy- 
namic form of electromagnetically induced transparency 
(EIT) [H, 1, % i, 0, S ©]. This light storage tech- 
nique has been shown experimentally to preserve opti- 
cal phase coherence [l(J as well as quantum correlations 
and states [ill . Il2| ; and thus has emerged as a promising 
technique for applications such as single-photon genera- 
tion on demand [H, 0, [H, [H, Il7l and quantum mem- 
ories [H, [H Hj and repeaters [H [H EH- However, 
practical applications will require significant improve- 
ments in the efficiency of writing, storing and retrieving 
an input photon state beyond values achieved to date 
[HI HH, HH, S3 ■ As an advance in this direction, we re- 
port in this Letter an experimental demonstration of an 
optimization protocol based on time reversal [28|, which 
determines the input signal pulse-shape that is written, 
stored, and retrieved with maximum efficiency for a given 
set of experimental conditions. This optimization proce- 
dure should be applicable to a wide range of ensemble 
systems in both classical and quantum regimes. 

We consider the interaction of a weak (classical or 
quantum) signal pulse and a strong (classical) control 
field with a A-typc medium. While the optimization pro- 
cedure is applicable to a wide class of regimes, we focus 
here on resonant light-atom interactions under EIT con- 
ditions. The group velocity of the signal pulse is pro- 
portional to the intensity of the control field, such that 
the quantum state of the signal pulse can be reversibly 
stored in a collective spin coherence (spin-wave) of the 
atomic ensemble by reducing the control field to zero [l| . 
For ideal writing, storage, and retrieval, the signal pulse's 
frequency components must fit well within the EIT spec- 
tral window (Aweit) to avoid incoherent absorptive loss: 
i.e., l/t s <C Albeit — V~dv g /L Q, where t s is the tem- 
poral length of the signal pulse, v g is the group velocity 
of the signal pulse inside the EIT medium at full con- 
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FIG. 1: (color online) Schematic of the iterative signal pulse 
optimization procedure described in the text, (a) An input 
signal pulse £i n (t) is mapped into a spin- wave S(z) using a 
control field envelope tt(t). (b) After a storage period r, the 
spin-wave is mapped into an output signal pulse £ out using 
the time-reversed control field envelope f2(r — i). (c) The 
time-reversed and normalized version of the measured £" ou t is 
used as the input £{ n in the next iteration. 



trol field intensity, L is the length of the medium, and d 
is the optical depth of the medium for the signal pulse 
in the absence of EIT conditions [2{|. In addition, v g 
must be small enough for the entire signal pulse to be 
spatially compressed into the ensemble before storage — 
i.e., v g t s <C L — so as to avoid "leakage" of the front edge 
of the pulse outside the medium before the back edge has 
entered. Simultaneous satisfaction of both these condi- 
tions is possible only at very large optical depth d, i.e., 
at high density and/or large sample size. However, oper- 
ation at very large d can degrade EIT performance and 
shorten the spin-wave coherence lifetime due to radiation 
trapping, competing nonlinear processes, etc. Therefore, 
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practical, high-efficiency light storage will likely be per- 
formed at moderately large d and require optimization of 
the input signal pulse-shape to minimize absorptive and 
leakage losses. 

Recently, a procedure to determine the optimal input 
signal pulse-shape for a given optical depth and control 
field was proposed 28]. This optimization procedure is 
based on successive time-reversal iterations and shown 
schematically in Fig. [TJ First, for a given input control 
field with Rabi frequency envelope ft(t), a trial input sig- 
nal pulse with envelope £- m (t) is mapped into a spin- wave 
S(z) inside the atomic ensemble (writing stage). (£j n (i) 
and the input control field are taken to be non-zero over 
the time-interval [— T, 0].) In general there will be some 
absorptive and leakage losses during this writing process. 
After a storage period r, an output control field Q(t — t) 

— i.e., the time-reversed version of the input control field 

— is used to map S(z) back into an output signal pulse 
£out(£)j which leaves the medium and is measured (re- 
trieval stage). The input signal pulse for the next itera- 
tion is then generated with a pulse-shape corresponding 
to a time-reversed version of the previous output signal 
pulse and an amplitude normalized to make the energy 
of the pulse equal to a fixed target value. These steps 
are then repeated iteratively, using the same input and 
output control fields, until the shape of the output signal 
pulse on a given iteration is identical to the time-reversed 
profile of its corresponding input signal pulse. The result- 
ing signal pulse-shape is predicted to provide the highest 
write/store/retrieve efficiency possible for a given optical 
depth and control field profile, and should be applicable 
to both quantum and weak classical signal pulses. 

In the experiment reported here, we tested this opti- 
mization procedure and confirmed its three primary pre- 
dictions: 

1. The write/store/retrieve efficiency (the ratio of en- 
ergies carried by the retrieved and input signal pulses) 
grows with each iteration until the input signal field con- 
verges to an optimal pulse-shape. See Fig. [51 

2. For a given control field profile and optical depth 
d, the optimization procedure converges to the same in- 
put signal pulse-shape and the same maximum efficiency, 
independent of the initial (trial) signal pulse-shape. See 
Fig.H 

3. For a given optical depth, different control field 
profiles result in different optimal signal pulse-shapes 
but yield the same maximum efficiency, provided spin- 
coherence decay during the writing and retrieval stages 
is small. See Fig. |4j 

We performed these experimental demonstrations us- 
ing a standard Rb vapor EIT setup, similar to that de- 
scribed in Ref. [3(J. A cylindrical 7.5 cm-long glass cell 
containing isotopically enriched 87 Rb and 40 Torr Ne 
buffer gas was mounted inside a three-layer magnetic 
shield to reduce stray magnetic fields. The Rb vapor cell 
was typically operated at a temperature ~ 60°C, corre- 
sponding to a Rb vapor density ~ 2.5 x 10 11 cm -3 and 
an optical depth d ~ 9.0. Optical fields near the Rb Dl 



transition (795 nm) were used for EIT and light stor- 
age. These fields were created by phase-modulating the 
output of an external-cavity diode laser using an electro- 
optical modulator (EOM) operating at the ground state 
hyperfine frequency of 87 Rb (6.8 GHz). The laser carrier 



frequency was tuned to the 5 2 S!/ 2 -F 
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transition and served as the control field during light stor- 
age; while the high-frequency modulation sideband, res- 
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onant with the 5 2 Si/ 2 -F = 1 
served as the signal field. The amplitudes of the con- 
trol and signal fields could be changed independently by 
simultaneously adjusting the EOM amplitude and the to- 
tal intensity in the laser beam using an acousto-optical 
modulator (AOM). Typical peak control field and sig- 
nal pulse powers ~ 5 mW and 0.1 mW, respectively. 
The laser beam was collimated to a Gaussian cylindri- 
cal beam of relatively large diameter (~ 7 mm) and then 
circularly polarized using a quarter- wave plate before en- 
tering the vapor cell. The Rb atom diffusion time out of 
the laser beam (~ 7 ms) was long enough to have neg- 
ligible effects 3l|. Small, remnant magnetic fields were 
the leading source of spin decohcrence, with typical spin- 
wave decay time constants ~ 2 ms. We used relatively 
short pulses and storage times, such that spin decohcr- 
ence had a negligible effect except for a modest reduction 
of the efficiency of the storage process. 

Fig. [5] shows an example implementation of the iter- 
ative optimization procedure, using a step-like control 
field and a trial input signal pulse with a Gaussian pro- 
file. Some portions of the first input pulse were incoher- 
ently absorbed or escaped the cell before the control field 
was turned off; but a fraction was successfully mapped 
into an atomic spin-wave, stored for 400 /is, and then re- 
trieved and detected. This retrieved signal pulse-shape 
was used to generate a time-reversed and normalized in- 
put signal pulse for the next iteration. After a few itera- 
tions, both the input and output signal pulses converged 
to fixed profiles, with the write/store/retrieve efficiency 
increasing with each iteration and reaching a maximum. 
In general, different trial input pulses all converged to 
the same optimal signal pulse-shape (e.g., see Fig. [3]). In 
addition, systematic variation of the signal pulse-shape 
uniformly yielded lower efficiencies than the pulse-shape 
given by the optimization procedure. 

We performed similar optimization experiments for a 
wide range of control field profiles. Some example re- 
sults are shown in Fig. [4] In general we found differ- 
ent optimized signal pulse-shapes for different control 
field profiles; however, the optimized write/store/retrieve 
efficiency was independent of the control field profile. 
This observation is consistent with the theoretical pre- 
diction [28[ that optimal light-storage efficiency does not 
depend on the control field, but only on the optical depth, 
provided spin decoherence and other loss mechanisms can 
be neglected during the writing and retrieval stages. 

To compare our experimental results with theoretical 
calculations, we approximated the 16-level structure of 
the 87 Rb Di line with a single three-level A-system and 
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FIG. 2: (color online) Top: Example data for the signal pulse optimization procedure, using a constant control field during 
writing and retrieval (timing indicated by dashed lines) and a 400 ^ts storage interval. Left: input signal pulses £i n , labeled 
by the iteration number and beginning with a trial Gaussian input pulse (iteration "0"). Center: signal pulse leakage for 
each iteration. Right: output signal pulse £ out for each iteration. Inset: Write/store/retrieve efficiency determined for each 
iteration from the measured input and output signal pulses: J E^dt/ J £? n dt. Bottom : theoretical calculation of the signal 
pulse optimization procedure, using the model described in the text and the experimental conditions of the measurements in 
the top panel. Inset: Calculated write/store/retrieve efficiency. 
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FIG. 3: (color online) Example data illustrating that the it- 
erative optimization procedure converges to the same signal 
pulse-shape and maximum efficiency independent of the ini- 
tial signal pulse-shape, for a given control field profile (shown 
in inset) and optical depth (d=9.0). 

modeled stored light dynamics as follows P, H, [28| : 

(dt +cd z )£{z,t) = igVNP(z,t), (1) 
d t P{z,t) = -~/P(z,t)+igVNS(z,t) + 

ifl(t- z/c)S(z,t), (2) 
d t S(z,t) = -i B S(z,t) + Xl(t-z/c)P(z,t)(3) 

Here £ is the slowly-varying envelope for the signal field, 
P is the transverse polarization of the atomic ensemble 
on the optical transition driven by the signal field, 5* is 
the spin-wave envelope, and f2 is the control field Rabi 



frequency envelope. Also, gy/N is the coupling strength 
between the atomic ensemble and the signal field, where g 
is the corresponding onc-photon Rabi frequency and N is 
the number of atoms along the laser beam that are avail- 
able to participate in light storage after optical pumping 
by the control field; 7 is the decoherence rate of P, due 
primarily to buffer gas collisions; and 7 S is the spin- wave 
decoherence rate. To calculate fJ, we used the dipole ma- 
trix clement of the \F = 2, m F = 1} -> \F' = 2, m F = 2) 
transition, which was the dominant control field transi- 
tion for our experimental conditions (32j. We approx- 
imated the laser beam profile as a uniform cylindrical 
beam with the same diameter and total power as used in 
the experiment. 

We solved Eqs. analytically by adiabatically elim- 
inating P, an excellent approximation for our experimen- 
tal conditions. We then calculated the results of the it- 
erative optimization procedure: i.e., the output signal 
field (both signal leakage and stored/retrieved pulses) as 
well as the generated input signal pulse for each succes- 
sive iteration. For these calculations we used values for 
g\fN \ 7, 7s, and the trial input signal pulse and control 
field appropriate for the particular experimental condi- 
tions. Example results of these calculations are shown 
in the bottom panels of Fig. [2] The calculated output 
signal pulse-shapes are qualitatively similar to the ex- 
perimental results and converge to an optimal input sig- 
nal pulse-shape within a few iteration steps. The calcu- 
lated efficiencies for the optimization procedure, shown 
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FIG. 4: (color online) (a)-(d) Examples of optimal input signal pulses determined experimentally by application of the iterative 
optimization procedure, for different control field profiles. In all cases, the initial (trial) signal pulse has the same Gaussian 
pulse-shape and amplitude as the trial pulse used in the data shown in Fig. [2] (e) Experimentally measured and (f ) calculated 
write/store/retrieve efficiencies as functions of the iteration number. 



in Figs. [D and BJf), are in reasonable agreement with 
experiment. We also confirmed that the effects of inho- 
mogencous Dopplcr broadening were small for the buffer 
gas pressure used in our experiments, by repeating the 
calculations in a more realistic approximation that in- 
cluded Doppler broadening of Rb atoms as well as veloc- 
ity changing collisions with buffer gas atoms. 

In conclusion, we experimentally demonstrated an it- 
erative optimization procedure, based on time-reversal, 
to find the input signal pulse-shape that maximizes the 
efficiency of light storage and retrieval. We confirmed 
the three primary predictions of the theory underlying 
the optimization procedure [28| : (i) efficiency grows with 
each iteration until the input signal field converges to its 
optimal pulse-shape; (ii) the result of the optimization 
procedure is independent of the initial (trial) signal pulse- 
shape; and (iii) the optimal efficiency does not depend 
on the control field temporal profile. We also performed 
theoretical calculations of the light storage process and 



the optimization procedure, and found good qualitative 
agreement with the experimental results, thus supporting 
the interpretation that optical depth is the key figure of 
merit for light storage efficiency. The optimization proce- 
dure should be applicable to both classical and quantum 
signal pulses and to a wide range of ensemble systems. As 
one example, since pulse-shape optimization with weak 
classical light pulses can be straightforwardly performed, 
such optimization could be used to determine the tempo- 
ral profile of input quantum fields, for which mode-shape 
generation and measurement are much more difficult to 
carry out. Also, pulse-shape optimization of the kind 
demonstrated here in atomic ensembles could be appli- 
cable to other systems, e.g., photonic crystals [33j]. 
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dation, and the Smithsonian Institution. 
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